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Abstract The stellar IMFhas been found to be essentially invariant. That means, 
any detailed observational scrutiny of well resolved populations of very 
young stars appears to show the IMF to be of more or less the same 
shape independent of the metallicity and density of the star-forming 
region. While some apparent differences are seen, the uncertainties in- 
herent to this game do not allow a firm conclusion to be made that the 
IMF varies systematically with conditions. The IMF integrated over 
entire galaxies, however, is another matter. Chemical and photometric 
properties of various galaxies do hint at galaxial IMFs being steeper 
than the stellar IMF, as is also deduced from direct star-count anal- 
ysis in the MW. These results are however sensitive to the modelling 
of stellar populations and to corrections for stellar evolution, and are 
thus also uncertain. However, by realising that galaxies are made from 
dissolving star clusters, star clusters being viewed as the fundamental 
building blocks of galaxies, the result is found that galaxial IMFs must 
be significantly steeper than the stellar IMF, because the former results 
from a folding of the latter with the star-cluster mass function. Further- 
more, this notion leads to the important insight that galaxial IMFs must 
vary with galaxy mass, and that the galaxial IMF is a strongly varying 
function of the star-formation history for galaxies that have assembled 
only a small mass in stars. Cosmological implications of this are that 
the number of SNII per low-mass star is significantly depressed and that 
chemical enrichment proceeds much slower in all types of galaxies, and 
particularly slowly in galaxies with a low average star-formation rate 
over what is expected for an invariant Salpeter IMF. Using an invariant 
Salpeter IMF also leads to wrong M/L ratios for galaxies. The detailed 
implications need to be studied in the future. 
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1. The shape of the stellar IMF 

The number of stars in the mass interval m, m + dm is dN = £(m) dm, 
where £(m) is the stellar initial mass function (IMF). The logarithmic 
slope is r(log 10 m) = dlog 10 £L(l°gio fn)/d\og w m such that log 10 £( m ) = 
log 10 k — alog 10 m for the power-law form, £(m) = km~ a , and a(m) = 
—dlog 10 £(m)/dlog w m = 1 — T(m). The logarithmic IMF is £L(m) = 
ra ln(10) £(m), which is useful when counting stars in logarithmic mass 
intervals. 

The stellar IMF is the distribution of stellar masses that results from a 
single star-formation burst. Because essentially all stars form in clusters 
(Lada & Lada 2003) we expect to measure the stellar IMF in young star 
clusters before they disperse to the field of a galaxy. 

The stellar IMF is one of the most fundamental distribution functions 
of astrophysics, and consequently a huge effort has been invested into 
constraining its shape since its first formulation in 1954 by Salpeter 
(1955) (at the Australian National University in Canberra) as a single 
power-law with a = 2.35 for 0.4 < m/M & < 10, based on an early 
analysis of star-counts in the solar neighbourhood. A further milestone 
in this remarkable scientific enterprise is given by Miller & Scalo (1979) 
who made a large effort in constraining the IMF to mass ranges outside 
the Salpeter limits and who deduced that the stellar IMF flattens below 
0.5 Mq. Scalo (1986) re-considered this problem in what today remains 
the most significant piece of work in this topic by studying all of the 
then available observational constraints on local star counts, the shape 
of the stellar luminosity function, the mass-luminosity relation, stellar 
evolution for early-type stars and the vertical structure of the Milky 
Way (MW), and suggested the IMF to turn-down below about 0.4 Mq, 
which had important implications in trying to understand the nature 
and occurrence of unseen matter in the disk of the MW. An important 
improvement in understanding the shape of the IMF for low-mass stars 
was contributed by Kroupa, Tout & Gilmore (1993 and two prior papers) 
who proposed that the strong maximum near My ~ 12 in the luminosity 
function of solar-neighbourhood stars stems from an inflection in the 
mass-luminosity relation, which in turn results from the association of 
H2 molecules and the onset of full convection below about 0.4 Mq. This 
revised the IMF to be essentially flat below 0.4 Mq, but the inclusion 
of corrections for unresolved multiple systems and detailed modelling of 
star-counts with Malmquist bias and galactic-disk structure solved the 
disagreement between local and deep star counts and thereby increased 
the slope of the IMF below 0.5 Mq. The resulting form of the two-part 
power-law IMF for late-type stars (c*i ps 1.3 : 0.08 ^ m/M ^ 0.5; Q2 ~ 
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2.3 : 0.5 < m/M Q < 1) has been verified by Reid, Gizis k Hawley (2002) 
using revised local star-counts that incorporate HiPPARCOS distances. 
For m ^ 1 Mq the IMF is proposed by Reid et al. to have 03 = 2.5 — 2.8, 
while Scalo (1986) derived 03 2.7, which is supported by the work of 
Yuan (1992). Written out Jo 



(7^-) a ° ,mi ow <m<m H , 

«m) = fc< ffi~Z, x- a ' mH - m<m °' (1) 

' fe) fe) ,m <m<m 1 , 



with exponents 



ao = +0.3±0.7 , 0.01 < m/M < 0.08, 

on = +1.3 ± 0.5 , 0.08 < m/M Q < 0.50, 

a 2 = +2.3±0.3 , 0.50 < m/M© < 1.00, 

a 3 = +2.7±0.7 , 1.00 <m/M , 



(2) 



is the KTG93 IMF after extension to the sub-stellar mass range (Kroupa 
2001; 2002). This IMF is the galaxial IMF as it is derived from galactic 
field stars, and, as we will see further below, differs from the stellar 
IMF. The multi-power-law description is used merely for convenience, 
as it allows us to leave the low-mass part of the IMF unchanged while 
experimenting with different slopes above, say, 1M Q . Other functional 
forms have been suggested (Miller & Scalo 1978; Larson 1998; Chabrier 
2001), but these bear the disadvantage that the whole functional form 
reacts to changes in the parameters. Such parametrisations are useful 
for studying possible changes of the stellar IMF with cosmological epoch. 

To constrain the stellar IMF we need detailed observations of young 
populations in star clusters and OB associations. The hope is that young 
clusters are not yet dynamically evolved so that the initial stellar popu- 
lation becomes evident. However, as dynamical modelling of young clus- 
ters containing O stars shows, significant dynamical evolution is already 
well established at an age of 1 Myr as a result of (i) the dense condi- 
tions before gas removal such that the binary-star population changes 
significantly away from its primordial properties, and (ii) rapid cluster 
expansion as a result of violent gas expulsion (Kroupa, Aarseth & Hurley 
2001). This would lead to a systematic biasing of the observed stellar 
IMF against low-mass members if the cluster was significantly mass- 
segregated prior to gas expulsion (Moraux, Kroupa & Bouvier 2004). 
In addition to these difficulties come obscuration by natal gas and dust 
and the large uncertainties in deriving ages and masses given inadequate 
theoretical models of stars that are still relaxing from their accretion his- 
tory, rotating rapidly, are variable and active, and have circum-stellar 
material. 
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The task of inferring the stellar IMF is thus terribly involved and 
uncertain. Nevertheless, substantial progress has been achieved. Given 
that there are many star clusters and OB associations, a large body of 
literature has been amassing over the decades, and Scalo (1998) compiled 
the then available IMF slope vs mass data. 

An updated form of this plot is shown in Fig. 1, and as can be seen, the 
stellar IMF essentially follows the above galaxial IMF below 1 Mq (af- 
ter correcting the observations for unresolved multiple systems), whilst 
having a Salpeter slope above 1M . Striking is that Large and Small 
Magellanic Cloud data (solid triangles) are not systematically different 
to the more metal rich MW data (solid circles). A systematic difference 
between dense clusters and sparse OB associations is also not evident 
(Massey 1998). Notable is that the scatter is large above 1 Mq but con- 
stant and that it can be fitted by a Gauss distribution of a values centred 
on the Salpeter index (Kroupa 2002). The scatter can be understood 
to result from statistical fluctuations and stellar-dynamical evolution of 
the clusters and the dynamically evolved OB associations (Elmegreen 
1999; Kroupa 2001). These observations on the a-plot thus point to a 
remarkable uniformity of the stellar IMF, which can thus be summarised 
by the multi-power-law form of eq. 2 but with 



Some apparently significant deviations from this standard or canonical 
IMF do occur though (Fig. 2), and it will be the aim of observers and 
modellers alike to try to understand why the two clusters of similar age, 
the Pleiades (Hambly et al. 1999) and M35 (Barrado y Navascues et al. 
2001), appear to have such different mass functions below 0.5 Mq. 

The existence of a universal canonical stellar IMF is therefore sup- 
ported by the majority of data. It has a Salpeter index above about 
0.5 M Q . But this result uncovers a possibly unsettling discrepancy be- 
tween this stellar IMF (eq. 3) and the galaxial IMF (eq. 2) which was 
found to be steeper for m > 1 Mq. 



The distribution of stellar masses in an entire galaxy results from the 
addition of all star-forming events ever to have occurred, 



= 2.3. 



(3) 



2. 



The galaxial IMF 




(4) 



This is the integrated galaxial IMF (IGIMF), i.e. the IMF integrated 
over space and time. Here £ ec i is the MF of embedded clusters, and 
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Figure 1. The alpha plot compiles 
measurements of the power-law index, 
a, as a function of the logarithmic stel- 
lar mass and so measures the shape of 
the MF. The canonical IMF (eq. 3) 
is represented by the thick short- 
dashed horizontal lines, and other 
functional forms are shown using other 
line-types (MS indicates the Miller 
& Scalo (1979) log-normal form; and 
the thin short- and long-dashed lines 
are from Larson (1998), while the 
thick short-dash-dotted curve is the 
Chabrier (2001) form labelled by Ch). 
Shaded regions indicate mass-ranges 
over which derivation of the MF is par- 
ticularly hard. For details see Kroupa 
(2002). 
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Figure 2. The measured logarith- 
mic stellar mass functions, £l, in the 
Orion Nebula Cluster [ONC, solid cir- 
cles], the Pleiades [triangles] and the 
cluster M35 [lower solid circles]. The 
average Galactic field single star-star 
IMF (eq. 2) is shown as the solid line 
with the associated uncertainty range. 
For details see Kroupa (2002). 



-^ecl,min = 5 M & , M ec i max ^ 10 6 Mq are the minimum and maximum 
cluster masses, respectively. Note that the IGIMF becomes indistin- 
guishable to the field-star IMF in galaxies in which presently on-going 
star-formation contributes insignificantly to the already present stellar 
population, and also that £igimf does not correspond to the shining 
matter distribution for m ^ 1 Mq . To evaluate this we need to consider 
only the recently formed stars. 
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Within each cluster stars are formed following the canonical IMF. 
However, small star-forming cloud cores, similar to the individual groups 
containing a dozen M & in gas seen in Taurus-Auriga, for example, can 
never from O stars. Since the stellar IMF has been found to be es- 
sentially invariant, even when comparing the small groups of dozens of 
pre-main sequence stars in Taurus-Auriga with the rich Orion Nebula 
Cluster (Kroupa et al. 2003), constructing a young cluster population 
is mathematically equivalent to randomly sampling the canonical IMF. 
However, sampling without a mass constraint would in principle allow 
small star-forming regions to form very massive stars, which would be in 
violation to the observations according to which massive stars are formed 
in rich clusters. We therefore choose stars randomly from the canonical 
IMF but impose a mass-limit on the mass in stars that can form from 
a molecular cloud core with mass M core = M ec \/e, where e^0.4 is the 
star- formation efficiency (Lada & Lada 2003), 

/■Umax 

M cc i = / m ■ £(m) dm, (5) 

where m\ ow = 0.01 M & . In each cluster there is one single most-massive 
star with mass m max , 

/•Tl max » 

1 = / £( m ) dm, (6) 

where m max * « 150 M & is the fundamental upper stellar mass limit 
(Weidner & Kroupa 2004a) . This set of two equations needs to be solved 
numerically to obtain m max = m max (M ec i) that enters eq. 4. We note 
that our usage of an m max (M cc i) correlation disagrees with Elmegreen's 
(2004) conjecture that no such relation exists. This problem is dealt 
with further by Weidner & Kroupa (2004b). 

The slope, (3, of the star cluster IMF, £ e ci oc M~j , is constrained by 
observations: 

P pa 2 locally (Lada & Lada 2003); 
(3 « 2.2 for LMC and SMC (Hunter et al. 

(3 ~ 2 ± 0.08 for Antennae clusters (Zhang 

Assuming (3 ~ 2.2 and the canonical IMF slope a = 03 = 2.35 for 
stars above 1 M we get the resulting IGIMF shown in Fig. 3 which 
is considerably steeper than the Salpeter IMF (see Kroupa & Weidner 
2003 for more details). 



■ 20 < M ecl /M < 1100: 

■ 10 3 < M eci /M Q < 10 4 : 
2003) 

■ 10 4 < M ccl /M < 10 6 : 
& Fall 1999). 
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Figure 3. The dotted line is the 
canonical stellar IMF, £(m), in loga- 
rithmic units and given by the stan- 
dard four-part power-law form (eq. 3). 
The dashed line is Cigimf(ti) for /3 = 
2.2. The IMFs are scaled to have 
the same number of objects in the 
mass interval 0.01 - 1.0 M©. Note 
the turn down near m max * = 150 M 
which comes from taking the fun- 
damental upper mass limit explicitly 
into account (Weidner & Kroupa 2004; 
2004a). Two lines with slopes ai ino = 
2.35 and an ne = 2.77 arc indicated. 





Figure 4. The IGIMF power-law in- 
dex Qigimf {m > 1 Mq) as a function 
of the star-cluster MF power-law index 
13 for a = 2.35,2.7, 3.2. 



In Fig. 4 the slope of the IGIMF, chigimf for m > IMq, is studied 
for various values of (3 and a. The effect is more pronounced the steeper 
(larger (3) the cluster mass function is. 

This has a profound impact on the evolution of galaxies as can be 
deduced from Fig. 5. In the upper panel is shown the number of white 
dwarfs per star relative to the Salpeter IMF. For (3 values above 2 it drops 
considerably, and in the lower panel, were the number of supernovae 
of type II (SNII) per star is plotted, the effect is even stronger. For 
a = 2.35 and (3 = 2.2 there are 89% of white dwarfs but only 35% of 
SNII compared to a Salpeter IMF. 



8 



0.95 - 




\\\ 

\\\ 
\ \ 
W 
\ 




SO. 85 - 



0.8 - 



\ N 

V 



0.75 - 



1 — 



0.8 - 




\ <>. 
\ o, 
Vv, 
\\\ 
\ \\ 
\ NX 
\ V- 



■B 0.6- 

CD 
L-. 



to 4 . 



0.2 - 







1.5 



2 



2.5 



Figure 5. Upper panel: The ratio of the number of stars in the IGIMF (for in- 
put stellar IMFs with a — 2.35, 2.7, 3.2, dotted, short-dashed and long-dashed lines, 
respectively) relative to the canonical IMF (Salpeter above 0.5 M©) in the mass in- 
terval 0.8 < m/M© < 8 (the relative number of white dwarfs) in dependence of 
the cluster mass-function power-law exponent (3. Lower panel: The same except for 
8 < m/Ai© < m max « = 150 M© (the relative number of SNII per star). The nor- 
malisation of the IMFs is as in Fig. 3. Note that the panels have different vertical 
scales. 

3. Variations among galaxies 

The star-formation rate (SFR) of a galaxy is given by 



ter IMF, £ ec i, is also invariant, and that M cc i jm i n = 5M & is fixed while 
M ec \ max may vary. Then the M ec i imax vs SFR relation can be calculated 
for different St. A comparison of this simple theory with empirical data 
is provided in Fig. 7. The figure plots empirical maximum star-cluster 
masses vs galaxial global SFRs, and a fit to these data yields 

log 10 (M eclmax ) = log 10 (fc M L) + (0.75(±0.03) ■ log 10 SFR) + 6.77(±0.02), (8) 

where /cml is the mass-to-light ratio. Fig. 6 shows this relation. Inter- 
estingly, the theory is nicely consistent with the data for (3 w 2.2 and 




Variation of the IMF 



9 



St fa 10 Myr, which suggests that the cluster IMF may indeed be quite 
invariant, and that a galaxy is able to assemble complete star-cluster sys- 
tems within typically 10 Myr independent of the SFR (Weidner, Kroupa 
& Larsen 2004). 

With the use of eq. 8 the time-dependent IGIMF of galaxies of dif- 
ferent types can be calculated: On specifying a SFR, M ec i max follows 
(eq. 8) and from eq. 4 the galaxial IMF can be computed in dependence 
of the star-formation history by adding up all galaxial IMFs generated 
in each multiple-<5t-epoch until the present. As a straight-forward result 
we expect the galaxial IMF to be steeper (larger chigimf) for low- mass 
galaxies than for massive galaxies because the average SFR is lower for 
the former. We also expect a dependence of the galaxial IMF on the 
star-formation history (SFH) of a galaxy. 

These expectations are born out. The final IGIMF of a dwarf or low- 
surface brightness galaxy with a stellar mass of M ga i = 10 M© assuming 
a stellar IMF slope a = 2.35 is shown in Fig. 7, and in Fig. 8 for a = 2.70. 
In the case with a = 2.35 three SFHs are considered. A single burst of 
star- formation followed by no further formation (solid line), an episodic 
SFR with 100 Myr long peaks every 900 Myr (long-dashed line) and a 
constant SFR over 14 Gyr (dotted line) . The influence of the SFH on the 
IGIMF is significant for galaxies with a small mass in stars. For a steeper 
input IMF (Fig. 8) this effect is even more pronounced. Such a steep 
IMF slope may possibly be the true value if the observationally derived 
stellar IMFs for massive stars are corrected for unresolved binaries (Sagar 
& Richtler 1991). In galaxies with a large stellar mass, M ga i = 1O 1O M 
(Fig. 9), the sensitivity of oigimf on the SFH becomes negligible because 
the average SFR is always large enough to sample the star-cluster IMF to 
the fundamental stellar upper mass limit. In summary, the IGIMF slope, 
ciigimf, is shown as a function of M ga i in Fig. 10. The differently shaded 
regions are for single burst, episodic and a constant SFR, respectively. 

What implications beyond a reduction of the number of SNII do these 
findings have? The observed diversity of metallicities in different dwarf 
galaxies of rather similar mass (Mateo 1998; Garnett 2004) may be ex- 
plained by this effect without invoking individually fitted effective yields 
as are needed in chemo-dynamical models with 'standard' IMFs (Lan- 
franchi & Matteucci 2004). On the other hand, our models indicate that 
massive spirals and ellipticals should show less variations in the IGIMF 
and thus metallicity due to their high average SFR which is actually 
what is observed in nearby galaxies (fig. 9 in Garnett 2004). Current 
available models for the chemical evolution of galaxies have difficulties 
to reproduce the chemical abundances in disc galaxies with a standard 
galaxial Salpeter IMF, as it produces too many metals (Portinari et al. 



10 



2004). This problem may be resolvable by the model presented here but 
detailed chemo-dynamical calculations incorporating our approach need 
to be computed. Finally, the mass-to-light ratia of galaxies are too low 
for models that assume invariant galaxial Salpeter IMFs; the real M/L 
ratia are larger. 




!) 



-6 -4 -2 

log 10 (SFR/[M s /yr]) 



8 I — i — i — i — i — i — r^i — i — i — i — i — i — i — i — | — i — i — i — i — r 




-2 - H 

-3 . I . . . . ' ■■■■■■■■ ] 

5 10 15 

t / Gyr 



-2-10 1 2 

log 10 (m/M s ) 



Figure 6. Maximum cluster mass 
versus global star-formation rate 
(SFR), both in logarithmic units. 
Filled dots are observations by Larsen 
(2001; 2002) with error estimates and 
the linear regression fit is the solid 
line (eq. 8). The other curves are 
theoretical relations (eq. 7) which 
assume the entire young-cluster 
population forms in 8t — 1, 10 and 
100 Myr (bottom to top). The cluster 
IMF has (3 = 2 (dotted curves) or 
(3 — 2.4 (dashed curves). 



Figure 7. IGIMFs for three dwarf or 
low-surface brightness galaxies with a 
final stellar mass of 10 7 M Q but differ- 
ent SFHs. The solid curve results from 
a single initial 100 Myr long burst of 
star formation. The dashed curve as- 
sumes a periodic SFH with 14 peaks 
each 100 Myr long and 900 Myr qui- 
escent periods in between. The thick 
dotted curve assumes a constant SFR 
over 14 Gyr. For all cases the canoni- 
cal power-law slope above 1 Mq (a — 
2.35) is used and the cluster IMF slope 
is (3 — 2.35. The straight solid line 
above 0.5 Mq shows the canonical in- 
put IMF for comparison. The thin 
dotted line is a Salpeter IMF extended 
to low masses. Note the downturn of 
the IGIMFs at high masses. It results 
from the inclusion of a limiting max- 
imum mass, m max < m max », into our 
formalism. 
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Figure 8. As Fig. 7 but only for 
two cases of the SFH: the burst case 
(solid line) and the low continuous case 
(thick dotted line), assuming a = 2.70. 



Figure 9. As Fig. 7 but for a galaxy 
with stellar mass 10 10 M Q and canon- 
ical a = 2.35. The IGIMF is only 
shown for a burst SFH (solid curve) 
and for a continuous SFH (dotted 
curve) . 



4. Conclusions 

1. IGIMFs are steeper than the stellar IMF. 

2. Therefore there are significantly fewer SNII per G-type star per 
galaxy. 

3. Chemical enrichment and M/L ratios of galaxies calculated with 
an invariant Salpeter IMF are wrong. 

3. The IGIMFs vary from galaxy to galaxy. 

4. Dwarf galaxies show the largest variations: 3.1 ^ qigimf ~ 3.6 for 
«3,truc = 2.35 (Salpeter) and = 2.35. 

5. The sensitivity of the IGIMF on the SFH increases with as- The 
true value of a^true for the stellar IMF may actually be significantly 
larger than 2.35 : 

6. Understanding the corrections for unresolved multiples on the deriva- 
tion of the high-mass part of the stellar IMF is of fundamental 
importance! 
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Figure 10. The resulting IGIMF slopes above 1M are shown in dependence of 
stellar galaxy mass for different models, as indicated. The lower bounds of the shaded 
areas are for an initial star formation burst forming the entire stellar galaxy while the 
upper bounds are derived for a continuous SFR. The various symbols correspond to 
calculated models. Symbols lying within the shades area correspond to models with 
an episodic SFH. 

Acknowledgement: 

This work is being supported by DFG grant KR1635/3-1. 
References 

Barrado y Navascues D., et al., 2001, ApJ, 546, 1006 
Chabrier G., 2001, ApJ, 554, 1274 
Elmegreen B.G., 1999, ApJ, 515, 323 
Elmegreen B.G., 2004, MNRAS, 342 

Garnett D.R., 2004, In: Cosmochemistry. The melting pot of the elements. XIII 
Canary Islands Winter School of Astrophysics, Tenerife, Spain, eds. C. Esteban, 
R.J. Garci Lopez, A. Herrero, & F. Sanchez, Cambridge University Press, 171 

Hambly N.C., et al. , 1999, MNRAS, 303, 835 

Hunter D. A., et al., 2003, AJ, 126, 1836 

Kroupa P., 2001, MNRAS, 322, 231 

Kroupa P., 2002, Science, 295, 82 



Variation of the IMF 



13 



Kroupa P., Bouvier J., Duchene G., Moraux E., 2003, MNRAS, 346, 354 

Kroupa P., Tout C. A., Gilmore G., 1993, MNRAS, 262, 545 

Kroupa P., Weidner C, 2003, ApJ, 598, 1076 

Kroupa P., Aarseth S., Hurley J., 2001, MNRAS, 321, 699 

Lada C. J., Lada E. A., 2003, ARAA, 41, 57 

Lanfranchi G.A., Matteucci F. 2004, MNRAS, in press (astro-ph/0403602) 
Larsen S.S., 2001, IAU Symposium Series, 207 
Larsen S.S., 2002, AJ, 124, 1393 
Larson R.B., 1998, MNRAS, 301, 569 

Massey P., 1998, in ASP Conf. Ser. Vol. 142, The Stellar Initial Mass Function, eds. 

G. Gilmore & D. Howell (San Francisco: ASP), p. 17 
Mateo M., 1998, ARA&A, 36, 435 
Miller G. E., Scalo J. M., 1979, ApJS, 41, 513 

Moraux E., Kroupa P., Bouvier J., 2004, A&A, in press (astro-ph/0406581) 

Portinari L., Sommer-Larsen J., Tantalo R. 2004, MNRAS, 347, 691 

Sagar R., Richtler T., 1991, A&A, 250, 324 

Salpeter E. E., 1955, ApJ, 121, 161 

Scalo J. M., 1986, Fundam. Cosmic Phys., 11, 1 

Scalo J.M., 1998, in ASP Conf. Ser. Vol. 142, The Stellar Initial Mass Function, eds. 

G. Gilmore & D. Howell (San Francisco: ASP), p. 201 
Reid I. N., Gizis J. E., Hawley S. L., 2002, AJ, 124, 2721 
Weidner, C, Kroupa, P., 2004, MNRAS, 348, 187 

Weidner C, Kroupa P., 2004a, Evidence for a fundamental stellar upper mass limit 

from clustered star formation, these proceedings 
Weidner C, Kroupa P., 2004b, Monte-Carlo experiments on star cluster induced 

integrated- IMF variation, these proceedings 
Weidner C, Kroupa P., Larsen S. S., 2004, MNRAS, 350, 1503 
Yuan J.W., 1992, A&A, 261, 105 
Zhang Q., Fall S. M., 1999, ApJ, 527, L81 



